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ABSTRACT 

This thesis develops s model for the estimation of 
the surface fluxes of momentum, heat, and moisture of the 
cloud topped marine atmospheric boundary layer by use of 
satellite remotely sensed parameters. The parameters 
chosen for the problem are the integrated liquid water 
content, q^, the integrated water vapor content, q y ^ , the 
cloud top temperature, the sea surface temperature, and 
either a measure of the 10 meter neutral wind speed or the 
friction velocity at the surface. 

Under the assumptions of a horizontally homogeneous, 
well -mixed boundary layer, the model calculates the 
equivalent potential temperature and total water profiles 
of the boundary layer along with the boundary layer height 
from inputs of q^, q yi « and cloud top temperature. These 
values, along with the 10m neutral wind speed or friction 
velocity and the sea surface temperature, are then used to 
estimate the surface fluxes using the methods of Stage 
(197?) and Liu, Katsaros, and Businger (1979). 

The development of a scheme to parameterize the 
integrated water vapor outside of the boundary layer for 
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the cam of cold air outbreak and California coaatal 
atratua ia preaantad. Tha achaina involvaa tha maan 
profilaa of ralativa humidity and tamparatura along with 
tha maan Jump of potantial tamparatura through tha 
invaraion to calculate tha integrated water vapor. In 
turn, this laada to a method of calculating the radiative 
tamparatura f lux at tha cloud top by uaa of Staley and 
Jurica'a (1970y 1972) radiation model. 

Sanaitivity atudiaa are preaantad ahowing tha 
potantial accuracy of tha technique. Tha model 'a expected 
accuraciaa with exiating aatellite technology are 24X far 
the heat flux and B2X for tha moiature flux. However, 
theae expected errora will be reduced with improvement of 
the aatellite parameter retrieval algorithma. 

In cancluaion, with improvementa of current 
technology, tha modal will provide raaaonabla eetimatea of 
the aurfaca fluxea. Tha modal ia tha only ona of ita kind 
currently available to eatimate tha aurfaca fluxaa of 
momentum, heat, and moiature of tha cloud topped marina 
atmoapharic boundary layer from aatellite moaaurabla 
paraincteri. 
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CHAPTER 1 i INTRODUCTION 


i . 1 PURPOSE 

The purpose of this thesis is to develop s model for 
the estimation of the surface fluxes of momentum, heat, 
and moisture by use of satellite remotely sensed 
parameters of the cloud-topped marine atmospheric boundary 
1 ayer . 

The thesis first introduces the problem. Second, the 
equations used to solve the problem are shown. Third, the 
computational form of the model is presented. Fourth, the 
results of the test runs of the model are given. Last, 
the thesis conclusions are drawn. 

The satellite parameters chosen for this problem are 
the Integrated liquid water content of the atmosphere, the 
integrated water vapor content of the atmosphere, the 
cloud top temperature, the sea surface temperature, and 
either a measurement of the 10 meter neutral wind speed or 
the friction velocity at the surface. This set of 
parameters is not the only possible set, but the 
technology exists today to remotely sense this set| hence. 
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it wai chosen for this study. The current technology to 
measure this set is discussed later in this chapter, but 
first the importance of this study is considered. 

1.2 SCIENTIFIC IMPORTANCE OF COLD-AIR OUTBREAK 

Modelling of the air-sea surface fluxes of momentum, 
heat, and moisture by use of remotely sensed satellite 
parameters is of importance due to the effect these fluxes 
have on the regional and global weather. 

Cold air outbreak has a pronounced effect on the 
western sections of mid latitude oceans off the east coasts 
of Asia and North America. In these regions, high winds, 
large air-sea temperature differences, and relatively dry 
air combine to give heat and moisture fluxes which are 
much greater than oceanic means. These fluxes over the 
warm water currents of the Kuroshio and Gulf Stream are 
known to have profound effects on cyclogenesis and storm 
intensif i cat ion, particularly on the development of 
wintertime cyclones from Taiwan and Cape Hatteras lows 
(Chou and Atlas, 1982| Agee and Howley, 1977). 

Recently, a method has been suggested by Chou and 
Atlas (1962) and Stage (19B3) to estimate the heat and 
moisture surface fluxes in the region between the shore 
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and the edge of a cloud bank during cold air outbreak, by 
remotely sensed parametara. This method uses ths boundary 
layer model developed by Stage and Businger (1961 a f b) for 
cold air outbreak. Further research is needed | however, 
to extend the estimation of the fluxes into the cloudy 
region. Hence, the reason for this study is the need of a 
model to estimate the surface fluxes in the cloudy region. 


1 . 3 THE MODEL 


The model developed in this study is capable of 
estimating the surface fluxes of momentum, heat, and 
moisture during a cold air outbreak by use of remotely 
sensed parameters of the marine atmospheric boundary 
layer. This model is based on the models of Stage (1979) 
and Liu, et aL (1979). It assumes unstable boundary 
layer conditions and uses the Businger diabetic profiles. 
Because of the assumption of unstable conditions, the 
model is not valid in neutral or stable conditions. 
However, the large air — sea temperature differences common 
in cold air outbreak or the large cloud top radiative 
temperature fluxes seen in California coastal stratus 


cases both produce the well -mixed boundary 
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condition* needed for the modal to run. The modal ia also 
limitad by tha praaanca of fog. Although tha modal ha* 
baan ravisad to includa tha praaanca of advaction fog, it 
ia atill unci aar aa to how tha aquation* which calculata 
tha aurfaca fluxa* raact in foggy condition*. Hence, fog 
ia an area in which furthwr atudy ia nacaaaary. 

1.4 SATELLITE TECHNOLOGY 

The technology ex lata today to meaaure tha 
parameter** cloud top temperature, integrated water vapor 
content, integrated liquid water content, aaa aurfaca 
temperature, and either a 10m neutral wind apead or a 
friction velocity for uae in thia model to aatimata tha 
aurfaca fluxea. 

Cloud top temperature may ba obtained from infrared 
aatallita measurements, such as the infrared sensors on 
the G0E8 east and west satellites. The other parameters 
may be calculated by uae of a scanning multichannel 
microwave radiometer (8HMR) such as the one carried aboard 
SEASAT or the one currently aboard NIMBUS 7. 

The 8MMR instrument 'a primary purpose ia to meaaure 
sea surface temperature and to provide a measurement of 
the surface wind speed.. However, SMMR also provides 
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estimates of atmospheric water vapor and liquid watar 
contanta. Tha 1 attar two paramatara ara found by acanning 
ovar varioua channala and employing an algorithm davalopad 
by Chang and Wilheit (1979). An survey of 8MMR*s 
capabilities may be found in the book, Microwave 
Ranjot® gana^QQ, edited by T.O. Allan (1983). 

As one can see, with existing technology the 
parameters needed for this model may be estimated from 
satellite measurements. In the next chapter, the model 


development is discussed. 
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CHAPTER 2 I MODEL DEVELOPMENT 
2.1 INTRODUCTION 

The modal development la discussed in thia chap tar. 
Aa atatad in tha firat chapter, tha purpoaa of tha modal 
ia to aatimata tha aurfaca fluxes of momentum, heat, and 
moiatura. Ona mathod of accompl lshing thia task ia to 
find an aatimata of tha variablaa 0 >f q^, and (Saa 

Figura 2.1.1). Thia may ba dona with aatallita 
maaauremanta of ©£, *^li» and °lvi * * nt * aasuming an 

unstable, wall mix ad boundary layar with and being 
conaarvad quantities <8ee Figura 2.1.2). Than, with an 

estimate of tha neutral wind spaed or tha friction 
velocity and tha saa surface temperature, the aurfaca 
fluxes may ba found by using tha bulk aerodynamic 

formulas (See Figura 2.1.3). 

Firat in this chaptar, aoma basic variablaa and 
thermodynamics are introduced. Second, tha derivation of 
aquations for tha modal are shown. Third, tha modal 
aquations are inverted. Fourth, tha aquations used to 
estimate tha aurfaca fluxaa are shown. Fifth, tha 
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parameterization of the liquid and Mater vapor contente 
outside of the boundary layer ie presented. Last, a 
method to calculate the radiative temperature flux at 
cloud top ie discussed. 


2.2 BASIC VARIABLES AND THERMODYNAMICS 

The basic variables used in this paper are non 
introduced. 

I'd DeMpoint temperature 

0 Potential temperature 

Sea surface temperature 
0^ Cloud top temperature 

Total water mi Mi no ratio 
q y ^ Water vapor mixing ratio at cloud top 
q^ Liquid water mixing ratio at cloud top 
z^ Cloud top height 

z £ Cloud base height 

z b Saturation level height 

k von Karmen's constant 

£ Ratio of molecular weights of water and air 

The basic thermodynamic properties are now shown. 
The equivalent potential temperature may be approximated 
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by 


© * <L v /C p’ q v 


( 2 . 2 . 1 ) 


whara L v la tha latant haat of avaporation for watar and 
C p la tha apaclfic haat of air at conatant praaaura. Tha 
total watar mi Ming ratio la dafinad aa 

q t -q y +q 1 (2.2.2) 

Tha Claualua-Clapayron aquation may ba writtan aa 

da-C (aL v £)/( (T d a )R) DdT d (2.2.3) 

whara a la tha watar vapor partial praaaura and £ la takan 

aa 0.622. Tha watar vapor partial praaaura la ralatad to 
tha mixing ratio by 

q y -(C/p)a (2.2.4) 

whara p^atmoapharic praaaura. Tha diffarantial of (2.2.4) 


gi vaa 
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(dq v /q v )-<d«/«)-(dp/p) (2.2.3) 

Substituting the hydrostatic aquation dp/dz“-^g, ths 
par fact gas law, p* p RT , and (2.2.3) into (2.2.3) givas 

(dq v /q y )-(L v e/( (T d a )R) )dT d +(g/(RT) )dz (2.2.6) 

Now, (2.2.6) may ba intagratad if ona knows tha 
saturation mining ratio, q r , corrasponding to soma 
dawpoint tamparatura, T dr , allowing tha dafining of two 
functional 

Q(T)-q r aKpC ( (-L y £)/R) (1/T-1/T dr ) 3 (2.2.7) 

and 

T daw <q> " C (R/<L v^ > >ln<q/q r >+1/T dr 3_l (2.2.8) 

Thus, Q(T) is tha saturation mining ratio for air at z»0, 
with tamparatura, T| and T d# ^(q) is the dawpoint for air 
at zrO, with water vapor mining ratio, q. These two 
functions are inverses! therefore. 
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■ T d.w (Q<T>> 


and 


<2. 2. 9a) 


q-Q(T . <q)> (2.2.9b) 

daw 

In ordar to avaluata thaaa functions. Stags (1979) 
obtainad saturation vapor prassuras from the Smithsonian 
Mataorological Tablas. Ha chosa valuas to giva tha bast 
fit for tha tamparatura ranga batwaan 0 and 10 “C. This 
givas tha function as 

Q (T) -5. 3542k 10- 3 {eKpC-5399. 286 ( 1 /T-3. 59505 k 10“*) 3> 

( 2 . 2 . 10 ) 

and 

T dBW (q)-5399.286/C21.0640-ln(qKl0 3 ) 3 (2.2.11) 

whara T and T dj#w ara in Kalvin and q and Q ara in g/kg. 
Tha liquid water potential tamparatura may ba dafinad as 

e i“ 6 a" (L v /C p >q t (2.2.12) 

and sines 0 m and q^ are conserved in non-pracipatating , no 
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ice, parcel motion, 0^ la also a conserved quantity. 
Below the saturation level, qj-O, and 0^ m 0 . The 
saturation level or lifting condensation level is 

expressed as 

z -(©.-T ..)/( T-(gT>/(L £>> (2.2.13) 

where A is the dry adiabatic lapse rate and (gT)/(L y £ ) 
may be expressed as the lapse rate of dewpoint. The last 
thermodynamic property to be introduced is the virtual 
potential temperature, 0 y . 

0 y - 0(l + (l/£ -l>q v - qi ) (2.2.14) 

Virtual potential temperature introduces the effects of 
both the buoyancy of water vapor and liquid water drag. 

2.3 DERIVATION OF MODEL'S EQUATIONS 

The equations used in this model were developed for 
the case of cold air outbreak. The marine atmospheric 
boundary layer (MABL) is assumed to be cloud topped and 
well mixed with 0^ and q^ being conserved quantities. The 
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cloud base i a aaaumad to ba aqual to the aaturatlon laval 
and that all mot Iona in tha cloud layer ara saturated. 
Tha cloud top is assumed to ba equivalent to tha base of 
tha inversion layar. 

A method is needed to determine values for the cloud 
top temperature, 6^, the integrated liquid water content, 
q li v * nc * the integrated water vapor content, q yi , of the 
marine boundary layer. The marine boundary layer model 
developed by 8tage (1979) provides such a method. Stage 
incorporates a routine to handle saturated partial 
behavior in which inputs of and z fa produce values of 
and q for a saturated parcel. Hence, with inputs of Q , 
q^. , and z^ it is possibls to find@ b , q^, and 

A method is now shown in which the cloud top 
temperature may be estimated by use of an iterative method 
with initial inputs of q^, and z^, and an initial 

guess for 6^. 

Cloud top temperature may be found by use of the 
following equations. 

T db“ e b" ( ^ -< 0 T>/<L vf >>2 b (2.3.1) 

The water vapor ml King ratio at cloud top may be found by 
using (2.2.9b) and (2.3.1) where 


PAGE IS 
QUALITY 
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q «Q(T ..) 
vb db 


(2.3.2) 


From (2.3.2) it follow* that 


dq vb -dQ/dTCd0 b -( * -(gT)/(l_ v C ) )d* b ] 


(2.3.3) 


Th* f lnit* difference form of (2.3.3) 1* 


dQ/dT«CQ(T db + & T)-q vb ]/AT 


(2.3.4) 


wh*r* A T i* a small increment. The quantity, C© bt i9 

expressed as 


^6 «C © -( ©.♦ 


(L../C ) q . 3 / C 1+ < L w /C« ) dQ/dT 1 


b + ' L v /c p ,q vb 


V- -p 


(2.3.3) 


Therefore, by making an initial guess at © b , one 
equations (2.3.1) to (2.3.3) and 


may use 


6 bn.w" 6 bold' , -^’ 6 b 


one may iterate numerically until a suitably small £&. i< 

b 

reached, thus providing © b and q vb . 
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The liquid water mi King ratio at cloud top may ba 
found from 

q lb“ q t -c, vb (2.3.6) 

In ordar to aolva for tha intagratad liquid water content, 
q lif tha depth of the cloud layer, d c must ba found. 
Using aquations (2.2.12), f or ©j , and (2.2.13), for z -f 
and assuming that i f , the saturation level, is equal to 
z c , tha cloud base level 

Ot-S-'c <2.3.7) 

Below cloud base, q^O. The liquid water mi King ratio is 
assumed to increase linearly with height from z £ to z fe , 
and than is assumed to ba equal to zero above tha boundary 
layer (Sea Figure 2.3.1). 

Thus tha q^ profile in tha cloud layer forms a 
triangle, and tha intagratad liquid watar content may ba 
found by 


q li -0.3Kq lb Kd c 


(2.3.6) 


ti ’ 


Tha intagratad total water content, q 


of tha boundary 



Sue fac e L*ve* 


Figure 2.3.1. Profile of in Marine Boundary Layer 




19 


layar i a axprassad as 



Sines is nail mix ad in tha boundary 1 ay ar 

q ti“ q t K1 b (2.3.9) 

Hanca tha final quantity, tha intagratad watar vapor 
contant, may ba found 

^vi"^ti~^li (2.3.10) 

Now that tha quantitlas, © fa , q Jif «nd q vi hava baan 
derivad, tha problam may ba invar tad. 

2.4 INVERSION OF EQUATIONS 


Tha pravious saction showad a mathod in which with 
inputs of 6^, q^, and war a usad to find tha satallita 
paramatars 0^, q ^ ^ , and • In ordar to computa tha 
boundary layar prof 11 as from thasa satallita paramatars a 
numarical mathod must ba found to invart tha sat of 
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aquations in ths last ssction. By invarting ths squat ions 
in Ssction 2.3 allows ons to taks ths satsllits msasurabls 
data and calculats z^ and ths boundary laysr prof i las of 
and q^. To bsgin, ths satsllits msasurabls paramatsrs, 
q li f q vi * * nd ars ussd as input. Than ons may 

calculats 


q ti“ q vi +q li 


(2.4.1) 


It is also possibls to sxprass q^ as 


•ti 


r *** r 

q t dz + 

Jt 


q t dz 


/■**? 

L -*• 


Taking / q^dz as a known quantity, this will ba 
discusssd in Ssction 2.6, and taking q^"*constant bslow z^ 
ons may writs 


'ti 


-f'. 


q t dz+q t xz b 


where 


r«i op 

q t xz b“ q ti -/ q t d * 

y *a 


(2.4.2) 


Ws now have ths quantity (q^xz^) from which q^ and z b may 
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be found by guaaalng and solving for q^. Than, by uaa 

of tha fallowing sat of aquations and a CNUlA fatal 
method, it ia possibla to estimate tha variables , q^, 

and which ara naadad in tha flux calculations. 

By guasaing z fa and solving 

q^“<q^xz b > /z fa (2.4.3) 

ona is laft with q^ and z fa . From aquation (2.3.1) one 

gats 

T db* 0 b“ < *-<0 T> '<«- v e >>* b (2.4.4) 

Tha water vapor mixing ratio may be found by 

t W Q<T db ) <2 - 4 - 5> 

Where Q(T) is the function described by aquation (2.2.10). 
Tha aqui valent potential temperature, is found from 


aquation (2.2.1) 



22 


ORIGINAL PAGE SS 
OF POOR QUALITY 


The liquid Mater mixing ratio at cloud top ia 


q lb" q t~ q vb 


<2.4.7) 


Tha dewpoint juat abova tha aaa aurfaca may ba found from 


T dO- T d.H <q t> 


<2.4.8) 


Whara T^^ <q > ia tha function axpraaaad by aquation 
<2.2.11). Tha liquid watar potential temperature cornea 
from aquation <2.2.12) 


W'vvt 


<2.4.9) 


Tha maturation level may be calculated uaing equation 
<2.2.13) 


*.-< 0 1 -T dO > /< »-<QT)/<L v ff )) 


<2.4.10) 


If ia graatar than z fa , no cloud ia formed in tha modal, 
and tha model anda. If z^ ia laaa than zero, than fog la 
praaant below cloud lavel. In tha caaa of fog, q^, ia no 


longer proportional to tha area of tha triangle ahown in 
Figure 2.3.1. Tha trapezoidal rule can now be applied aa 
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can be seen in Figure 2.4.1. 

The valua of at tha cloud base is cal lad Qj c > 
Under no fog conditions, qj c «0. In tha case of fog, tha 
cloud base ia at tha surface and q^ c is not equal to zero. 
Note that the term "fog" refers to advection fog. If fog 
is present in the boundary layer, q^ c may be found from 
the input values of 0^, q^., and z s by the same iterative 

method shown in Section 2.3. 

The dewpoint temperature at z^ may be found by 

T d> - 3-<gT>/<L v £ ) >z a (2.4.11) 

where a value for & a has been guessed. The water vapor 
mixing ratio at saturation level is 


q «Q(T. > 
^vs ds 


(2.4. 12) 


where Q(T> is from equation (2.2.10). A small increment, 
AT, is now introduced in which 


dsm 


■T . , .♦ 

iw dsold 


AT 


dsnew 


(2.4. 13) 


The mixing ratio of T 


is 
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■’ <T d.n M , ‘ 0,T d.n«,' <2 -‘ , - ,4 > 

found by tha maw mathod as (2.2.11). Thw finltw 
diffwrwncw form, aquation (2.3.4), ia now amployad 


dQ/dT.(q(T d-n<tw )- qv- )/ AT 


(2.4. 15) 


From aquation (2.3.3) 


fa 


c ©.-( © B +(L v /C p )q VB )3/Cl+(L v /C p )dQ/dT3 


(2.4. 16) 


If °& B i» la*a than aoma tolaranca laval than qj c may ba 
found by 


q lc" q t" q vs 


(2.4. 17) 


Otharwisa, aquations (2.4.11) to (2.4.16) ara rapaatad 
with 


G 

•n«M 


e 


ol d 



(2.4. 18) 


unti 1 


q lc is found. 

If 0<z <z. than a cloud is for mad in tha modal and no 
C D 
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fog is present. Then, the triangle rule illustrated by 
Figure 2.3.1 applies. The equation used to handle both 
fog and no fog cases is 

q li -0.3d c (q lb +q lc > <2.4.19) 

The difference between the value of q^ calculated in 
<2.4.19) and the value of q^ measured by the satellite is 
then checked to see if it is below an acceptable tolerance 
value. If it is not, then a cegyla falsi. routine is used 
to calculate a new value of and equations <2.4.3) to 

<2.4.19) are repeated until an acceptable value is 
achieved. Once this occurs, the boundary layer profiles 
of 0^ and q^ are found. It is then possible to begin the 
calculations to derive the surface fluxes of momentum, 
heat, and moisture, which is the topic of the next 
section. 

2.5 CALCULATION OF SURFACE FLUXES 

The calculation of the surface fluxes of momentum, 
heat, and moisture may now be done. As discussed in the 
first chapter, satellite measurements of both a neutral 
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Mind ipHd at some height and a friction valoclty may ba 
calculated by maana of ampirical formulas. Because the 
debate continues on Mhich is a more appropriate 
measurement in its correspondence to the measured 
backscatter, I have chosen to alloM the use of either in 
the model. With the addition of satellite measured sea 
surface temperature, the surface fluxes may be found. 

By folloMing the method of Liu, ft <1979), the 
surface fluxes may be estimated by first finding u # and 
Zq. Then, the roughness Reynolds number is computed from 
Mhich z^ and z^ are found. With the addition of the sea 
surface temperature and mixing ratio, the surface fluxes 
may be found by using the bulk aerodynamic formulas and 
the Busi nger-Dyer unstable diabetic profiles. 

For input involving a neutral stability Mind, Kondo's 
(1975) relation is used 

C dn“ p '*' q<u 10n > " <2.5.1) 

Where p, q, and r are constants Mhich vary for Mind speed 
velocities betMeen 0.3 and 50.0 m/s. Kondo's relation Mas 
chosen because of its ability to cover large ranges of 
Mind speed. This sIIoms the model to incorporate cases 
such as California coastal stratus Mhich do not exhibit 
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tha larga alr-ata tamparatura diffarancaa common for cold 
air outbreak p thus, craating a naar nautral atmoapharic 
atability. 

Equation (2.3.1) may ba axpraaaad aa 

< u # / u l0n )*-P^qC <u # /k)ln<10/* 0 ) 3- (2.3.2) 

whara I hava uaad 

C dn- (u * /u 10n >3 (2.3.3a) 


u 10n" <u * /k) <ln<1 ° / *0 >> (2.3.3b) 

Ona may calculata tha friction valocity, u # , by 

u .-< c dn (u 10n* ,> *~ <2.3.4> 


Tha roughnaaa langth, z^ f may ba found from tha axpraaaion 
for u }Q n in aquation (2.3.3b). _ 

* 0 -10/Caxp((u l0n k)/u*>3 (2.5.5) 


For input involving aatallita darivad friction 
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velocity, Charnock's (1959) relation ie used to find the 
roughness length. The equation is 

* 0 -a(u # ) a /g (2.5.6) 

where g is the gravitational acceleration and a is 
Charnock's constant. Wu (1969) in his study of sea 
surface roughness, found a*0.0156 to be a reasonable value 
for a large range of oceanic data. 

Following the method of Liu, e£ al & (1979) the 
roughness Reynolds number is expressed as 

R r -*0U # /9 (2.5.7) 

where V is the kinematic viscosity. 

According to Liu, et (1979), the variables z^ and 

Zq may be expressed in terms of the roughness Reynolds 
number by 

z t -(ajR r -‘ V >/u # (2.5.6) 


and 
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Z q" ( *2 R r M V >/u * (2.5.9) 

wh era the coefficients , bj , a^ and b 2 * or different 
ranges of R r are mhown in Table 2.5.1, taken from Liu, 
ali (1979). 

Tha distribution of valocity, tamparatura, and 
humidity for tha aurfaca iayar, outaida of tha ragion 
where molacular af facta dominata, are govarnad by tha 
Buaingar prof i las 

(u- u 0 )/u*-Cln<*/i 0 >- ^u 3/k (2.5. l0> 

( 0 - 0 o )/© # -Cln(z/z t )- ^ t 3/«< h k (2.5.11) 

(q-q 0 0/q # -Cln(z/z q )- ^3/^^k (2.5.12) 

Tha diabatic tarma, ^ u , / ¥' t , and ara givan by 

^ u -2 lnfCl-Kl-/ff u (z 0 /L) ) »~3/2>+lntCl + (l- y0 u (z o /L) ) *'*3/2> 
-2tan -1 C <1- /$ u (z 0 /L) ) *'*]+ V/2 (2.5. 13) 

1* t -2 ln<Cl + (l-^<z 0 /L> ) *'=*3/23 


(2.5. 14) 
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Tab la 2.9.1. Loaar Boundary Valuaa of tha Logarithmic 

Prof i la. z . u_ / -a f R^ b ‘ and z u m / -a-R 
t * ir q * £ r 


R r 


•l 


*2 

b 2 

0. - 0 

.11 

0. 177 

0. 

0.292 

0. 

0. 11 - 

0.825 

1 . 376 

0.929 

1.808 

0.826 

0.825 


1.026 

-0.599 

1.393 

-0.528 

3.0 - 

10.0 

1.625 

-1.018 

1.956 

-0.870 

10.0 - 

30.0 

4.661 

-1.475 

4.994 

-1.297 

30.0 - 

100.0 

34.904 

-2.067 

30.790 

-1.845 
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•nd 

^-2 ln«l-Ml-/> q (z q /L) > l '»J/2> (2.3.15) 

raapacti valy. Liu, *jfc eli (1979) uaad k«0. 4, 

l/(o< h k)-l/(«< w k>-2.2, and 

To includa tha affacta of tha buoyancy of watar 
vapor, tha Obukhov langth, L, nay ba axpraaaad in tarna of 
tha virtual potantial tamparatura flux 

L-< e v u # *)/(gk © v# ) 

wh « r * ® v coffin from aquation (2.2.14) 

& v - (l+(l/£ -l>q t > 

and 0 y# in axpraaaad aa 

© y# - 0 # (l-(l/£ -l)q t )^(l/e -1)0 q # 

Tha profilaa (2.3.10) to (2.5.12) ara formulatad in 
tarma of naan quantltiaa for tha boundary layar. Thaaa 
naan valuaa can ba aaaoclatad with tha aurfaca flux 
profilaa. Naar tha aurfaca, machanical ahaar domlnataa 


(2.3.16) 


(2.5.17) 


(2.5. 16) 
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and large gradients may exist in the mean profiles. 
Higher up, buoyancy begins to dominate and the shears 
gradually vanish as the profiles approach their mixed 
layer mean values. 

If a matching height, h, is Introduced it is possible 
to solve the problem of matching the mean mixed layer 
profiles with the surface layer profiles (Wyngaard, ft 
§L , 1974). Let h be the height at which the profiles 

(2.3. 10) -(2. 3. 12) produce the mixed layer mean values. 
Then one may write 

In (h/* 0 ) — (h/L) "In <-h/L> - / t* u <-h/L> -In <-i 0 /L) 

ln(h/* t >- (h/L)-ln(-h/L)- ^ (-h/L)-ln(-* t /L) (2.3.19) 

In (h/* )- tL (h/L) -In (-h/L) - V (-h/L) -In (-z„/L) 

q q q q 

The height at which the gradients become vanishingly 
small is determined by the relative importance of shear 
and buoyant production of turbulence kinetic energy and 
thus h should be proportional to -L (Wyngaard, flkiLf 

Therefore, 


1974) 
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K u -ln(-h/L>- ^ u <-h/L> 

K t -ln(-h/L>- \ <-h/L) (2.3.20) 

K q -ln(-h/L>- ^(-h/L) 

ara constants and ths prof lias (2.3. 10) -(2. 3. 12) may ba 

•xprassad as 

(u-u 0 >/u # -CK u -ln(-* 0 /L) 3/k (2.3.21) 

(©- © Q )/ © # -CK t -ln(-* t /L)3/o<. h k (2.3.22) 

(q-q 0 )/q # -CK q -ln(-z q /L) 1/ci ^k (2.3.23) 

Wyngaard, at al . (1974) took -h/L«-10.0 f thus giving 

K u — 0.2 and K t -K —1.02. 

In ordar to solva aquations (2. 3. 21) -(2. 3. 23) tha 
Obukhov langth, L v must ba found. Equations (2.3.16)- 
(2.5.18) show that it is nacassary to know 0 # and q # in 
ordar to do this problam. If an initial guass of © v# is 
givan and aquation (2.3.17) is usad to solva for d v , tha 
Obukhov langth may ba datarminad by aquation (2.5.16). 
Tha paramatar 6 # may ba found from aquation (2.3.22) and 
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the satellite measured surface temperature 6>^. 

© # -(0j- © 0 >/C * h k<K t -ln<-i t /L> ) 3 (2.5.24) 

The sea surface value of the water vapor ml King ratio* q Q , 
may be found by using the function Q(T) expressed in 
equation (2.2.10). 

q 0 «Q(© 0 ) (2.5.25) 

Assuming no liquid water or ice is present in the mixed 
layer where the mean profiles exist* q^ may be used as the 
value of the water vapor mixing ratio of the mixed layer. 
Thus, q # may be found from equation (2.5.23). 

q # -(q t -q Q )/C «* B k(K q -ln(-* q /L) )3 (2.5.26) 

A nw v.lu. of e y# may now be calculated using equation 
(2.5.18). The difference between the guessed value and 
the calculated value of 6 v# may be found 

A© - © - ft 

v* vtgueis v*calc 


(2.5.27) 
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Whan 00 reaches a small valus approx imataly aqual to 
zaro, © v# haa baan found, and therefore, tha values of © # 
and q # an known. 

Tha surface fluxas assoc 1 at ad with tha prof lias 
expressed in aquations (2. 3. 21 >-(2.3. 23) may ba writ tan as 

(2.3.28) 


(2.3.29) 

"'Q t '| 0 — u # q # (2.3.30) 

where u'w'l^ is tha surfaca flux of momentum, w' $ 'I q is 
tha surfaca flux of sanaibla heat, and w'q^'Jg 11 tha 
surfaca flux of molstura. Sinca u #f © and q # art now 
known, tha fluxas may ba calculated from thasa aquations. 

Tha method listed above for determining the fluxas 
assumes unstable conditions. Therefore, a limit must ba 
place on tha calculations for whan conditions approach 
near neutral stability. 

As conditions approach near neutral stability, tha 
relation z^/L goes to zero sinca L goes to plus or minus 
infinity. Tha author chose to limit z q/L to ba 


uw'lo— u*» 
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<* 0 /L> min £— l.OxlO-^ <2.5.31) 

By taking a typical valua of z^ ovar tha ocaan ai 2.0xl0 - ‘* 
matara (Mlyaka, at ai it 1970) v thia givaa an uppar limit 
of tha Obukhov langth of 

~L"2000 matara 

Tha maNimum haight of z b in tha modal la 2000 matara. In 
ordar for tha Obukhov langth not to axcaad z^, tha 
aaaumptlon axprassad in aquation (2.5.31) appaara to ba a 
valid lowar limit for naar nautral atability conditions. 

2.6 PARAMETERIZING WATER VAPOR ABOVE THE MABL 

A mat hod to paramatariza tha watar vapor contant 
abova tha MABL ia now praaantad. An analysis la naadad in 
which tha prof 11a of q^ abova tha boundary layar may ba 
found. From this prof i la, tha lntagratad watar vapor 
contant may ba found alnca qj is aasumad to ba zaro abova 
z^. A climatological profila ia tha aasiaat maana 
aval labia of accomplishing tha analyaia. Schamas for both 
cold air outbraak and California coastal stratus ara givan 
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in thla •action. 

2. 6. a Cold Air Outbreak 

The profile for cold air outbreak was derived from 
the NASA Hesoscale Air — Saa EXchanga (HA8EX) project. In 
ordar to find tha prof i la of abova tha boundary layer, 
prof i las of temperature, ralativa humidity, and mi Ming 
ratio wara studiad. Tha profiles of tamparatura and 
ralativa humidity wara chosan to find tha profile of q^.. 
Tha mining ratio profile was discarded because tha 
ralativa humidity profile seamed more readily adaptable to 
a wide range of situations. 

Given tamparatura and ralativa humidity, one may find 
So that tha profiles studiad in MASEX may be made 
relevant to other cases, tha lapse rates of tamparatura 
and ralativa humidity Mara found. In addition, the mean 
jump of tamparatura from tha bottom of tha inversion to 
tha top of tha inversion was calculated along with tha 
mean ralativa humidity at tha inversion top. By using 

these mean values of ^and RH , and tha lapse rates of 
0 and RH, tha profile abova tha boundary lay or of q^ may 
ba found. 

Tha mean values of and RH ara 3 "C and 4SX, 
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res pectively. The laps* rates of RH and 0 ara 8 X/km and 
6 “C/km, raapactivaly. By adding to tha cloud top 
temparatura and intagrating ovar 10m ataps using the lapse 
rates listed above, q yi above tha MABL may be found. 

2.6.b California Coastal 8tratus 

Similarly, a profile for California coastal stratus 
was computed. Analysis of tha profiles found in the 
studies of Naiburgar, gt al*. (1961) and lieitin (197S) 
yielded tha mean values for this case. Values ofA^«7“C, 
RH«35X, lapse rate of temper at ure*6 “C/km, and lapse rate 
of RH“B X/km were extracted from the data. 

2.7 RADIATIVE FLUX 

Given the profiles above the boundary layer, the 
radiative flux at cloud top may be found by using the flux 
emissivlty method of Staley and Juries <1970|1972>. Since 
the model assumes no cloud above the boundary layer, the 
flux at the top of the atmosphere is equal to zero. 

The optical path length, u(z), is given by the 
integration of the water vapor mixing ratio profile. 
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u(z> 



(q/g)dp 


<2.7.1) 


In equation (2.7.1), *top the highlit sounding 
lsvsl svsilsbls. Ths ■flux emissivity <u) is found by 
ths intsrpolstion of ths vsluss tsbulstsd by Stslsy and 
Juries} thus, giving ths downward radiativs tsspsraturs 
flux as 


(2.7.2) 

sines ths atmosphere is assumsd to bs cloud frss abovs ths 
MABL. 

Ths method dsscribsd abovs, or any othsr similar 
radiation modal , will ylsld a valua for ths radiativs 
tsspsraturs flux. Dus to ths uncertainties involved whan 
using mean profiles in place of actual profiles, the 
radiativs temperature fluxes calculated by the model are 
not yet reliable. When a better method of finding the 
profile of water vapor above the MABL is found, the 
radiative flux method listed above will be useful. 
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CHAPTER 3 « COMPUTATIONAL FORM OF THE MODEL 

3.1 INTRODUCTION 

Tha computational form of tha modal la diacuaaad in 
this chap tar. Tha aquations for tha modal ara givan in 

Chaptar 2| hanca, this chaptar rafarancas tha aquations 
instaad of rapaating tham. This chaptar shows tha staps 
naadad to produca astimatas of tha surfaca flux as from tha 
input data. 

3.2 THE MODEL 

In ordar to astimata tha surfaca fluxas of momantum, 
haat f and moistura, tha satallita inputs of cloud top 
tamparatura, intagratad liquid watar contant, and 
intagratad watar vapor contant must ba convartad into 
and tha boundary layar profilas of Q and q^. First a 
guass of tha boundary layar haight is givan. Than, tha 
intagratad watar vapor contant outsida of tha boundary 
layar is found by tha mathod in Saction 2.6. Than, tha 
mathod givan in Saction 2.4 is usad to calculata £>^, , 
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and z fa . The integrated liquid water ia calculated from 
tha values found in Section 2.4. This calculated value is 
then compared to the satellite measured value. If the 
difference bet w een the measured and calculated values of 
is approximately equal to zero, the model continues. 
Otherwise, a regula fals^. routine is used to obtain a new 
value -of the boundary layer height, and the process above 
is repeated until the difference between the measured and 
calculated q^'a Is approximately equal to zero. 

Upon finding the values of 0 -f q fcf «nd z b , the model 
continues. If a neutral wind speed is provided, the model 
calculates the friction velocity and roughness length by 
the method shown in equations (2.5. 1) to (2.5.5). On the 
other hand, if a friction velocity is given, the roughness 
length is calculated by equation (2.5.6). With the values 
of u # and Zq, and following the method of Liu, §i.i_ 

(1979) , a roughness Reynolds number may be computed. This 
is accomplished by equations (2.5.7) to (2.5.9). 

The last satellite measured value, the sea surface 
temperature, is now used. Values of 6 # and q # must be 
found so that the fluxes may be computed. From 0^ and q^. 


is found. The 

quantity (dj-dg) 

used 

for 

the 

computati oh 

°f e # . 

The quantity (q^-q^) 

, where q Q 

is 

found from 

the sea 

surface temperature. 

is 

used 

for 



43 


calculating q # . With an initial guess of 0 y# , aquations 
(2.3.16) to (2.5.18) and aquations (2.3.24) to (2.3.25) 
ara amploysd to calculata 0 # and q # . Than 0 y# is computad 
by aquation (2.5.26) and compared to tha valua of tha 
original guass of © y# by aquation (2.5.27). If tha 
diffaranca is not approximately aqual to zero, a re pu la 
falsi routina is employed to compute a nax 0 y# . Tha 
procass abova is rapaatad until tha diffaranca is 
approximately aqual to zero, at which time the modal 
continues. 

Finally, with values for u # , 0 # , and q # , tha surface 
fluxes of momemtum, heat and moisture may be estimated. 
Equations (2.5.28) to (2.5.30) ara used to accomplish this 
task. 

A note about tha computations listed abova. Original 
test runs of tha modal showed that more than one solution 
existed for each combination of input parameters. 
Numerically, these solutions ware valid| however, these 
solutions ware physically unrealistic due to tha extremely 
high dewpoint temperatures and low mixing ratios 
encountered. In order to “force" the model to converge to 
the “real istic" answer the following was done. An upper 
limit of 2000m for z fa was chosen, since the majority of 
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the physically unrealistic solutions were at z^ values 
greater than 2000m. To handle the unrealistic solutions 
that occurred below 2000m, the initial guesses of z^ for 
the regula falsi routine were set at 1200 and 100m. Since 
no double solutions were seen to exist below 1200m f any 
solutions with z^ less than this height will be accurate. 
For solutions between 1200 and 2000m v where double 
solutions might occur, a very restrictive tolerance was 
placed on q^. Runs of the model showed that was 
extremely sensitive to changes in height. By requiring a 
near perfect convergence of the measured and calculated 
integrated liquid water contents, the problem of a 
possible double solution in the 1200 to 2000m range seems 
to have been solved. 

In this chapter, the computational form of the model 
has been discussed. The steps needed to estimate the 
surface fluxes from the original satellite data have been 
shown. The next chapter deals with the model 's results 
for its test runs. 
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CHAPTER 4 i MODEL RESULTS 
4. 1 INTRODUCTION 

Tha rnulti of tha tost run* of tha modal an 
diicuaMd in this chaptar. First, ths applications of tha 
modal ara prasantad. Naxt, tha parcant variation of tha 
output paramatars with variation of tha input paramatars 
is shown. Naxt, tha ax pact ad accuracy of tha modal with 
raspact to currant accuracias of aatallita maasuramant* is 
prasantad. Last, tha strangths and waaknassas of tha 
modal ara discussad. 

4.2 MODEL APPLICATIONS 

Tha modal was originally davalopad for tha casa of 
cold air outbraaki howavar, it may also ba appliad to 
California coastal stratus. Tha modal assumptions of a 
horizontally homoganaous, cloud— toppad boundary lay ar with 
littla or no liquid watar abova tha boundary layar apply 
to both casa* mantionad abova. 

In ordar to tast tha modal 's capabilitias, tha 
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routina described in Saction 2.3 mu run ovar a larga 
ranga of valuaa for 0^, q^, and z^. Tha pradictad valuaa 

of 0 b f q^, and q yi from this routina war a than uaad aa 

input for tha modal. For tha caaaa in which tha routina 

in Saction 2.3 pradictad clouda and/or clouds and fog, tha 
modal computad valuaa of ©^, q^ , and z fa which agraad with 
tha variablaa inputtad into tha routina. Howavar, for 
caaaa in which no clouds wars pradictad by tha routina, 
tha modal would not con varga to an anawar. Tha routina in 
Saction 2.3 was than run for typical rangas of valuaa for 
tha cut of cold air outbraak with ranging from -30 to 
10 °Ci 3^ ranging from 1.0 to 3.0 g/kg, and z b ranging 
from 500 to 2000 matars. Tha routinas pradictad valuaa of 

0b* 3li* and q V £ wars than uaad aa tha ranga of valuaa to 

computa tha sansiti vitiaa of tha output paramatars 

deacribad in tha naxt saction. 

4.3 VARIATION OF INPUT AND OUTPUT PARAMETERS 

Tha parcant variation of tha haat and moiatura fluxaa 

with variation of u 1<)n , u # , ® 0 . q li* * nd q vi i- now 

praaantad. Tha modal was taatad ovar tha ranga of valuaa 
liatad abova plus tha rangaa for u 1Qn of 5 to 20 m/a, for 
u # of .1 to .7 m/a, and for of 5 to 20 “C. 


Tha 



47 


ORIGINAL PAGE IS 
OF POOR QUALITY 


flguras prnvnttd in this chapttr ara from ona tost rang* 

of valuas for tha modal. Tha input variables for tha 

cmmmu shown hara wars hald constant, sKcapt whan they wara 

tha variabla baing t as tad, at u 10n "10 m/s, 9 b — 13 ”C, 

© 0 "10 "C, qjj*300 m g/kg, and q y j-900 m g/kg. 

Tha 10« nautral wind spaad was tastad ovar tha ranga 

of 3 to 14 m/s in this casa. U JOn is ralatad to tha 

surfaca flux as through u # , z Q , z t , and z q . Tha rasults of 

tha variation of both tha haat and moistura flux as with 

raspact to is naarly linaar. Analysis shows that 

with a +/-10X variation in u 10n rasulting in a +/-10X 

variation in both tha haat and moistura fluxas. 

Tha friction valocity was tastad naxt with a ranga of 

.2 to .63 m/s in this casa. Tha friction valocity is 

dlractly ralatad to tha fluxas as wall as through z^, z^, 

and z . Tha rasults of tha tasts showad tha variation is, 
d 

as in tha casa of u^, naarly linaar. Analysis of tha 
rasults shows that a +/-10X variation in u # producas a 
+/-9X variation in both tha haat and moistura fluxas. 

Naxt, tha cloud top tamparatura was tastad ovar a 
ranga of -13 to -10.5 °C in this casa. Tha cloud top 
tamparatura is ralatad to tha haat and moistura fluxas 
through 0j and q^. Rasults for tha variation of tha haat 
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f 1 uk may bi seen in Figure 4.3.1. 

As Been in Figure 4.3.1, the heat fluK decreases fcr 
increasing 0 ^ . The change in slope of the curve is due 
to the presence of fog predicted by the model. Because of 
the uncertainty of the model *s evaluation of fluxes under 
foggy conditions, only cases in which no fog is present 
were considered for error determinations. The error for 
the no fog portion of this case as well as in all other 
cases is for a +/-1 "C change in 0 fa , the heat fluK varies 
by -/+4X. Figure 4.3.2 shows the variation of the 
moisture flux with 0^- 

As can be seen in the figure, the moisture flux also 
decreases with increasing The change in slope of the 
curve seen in the figure corresponds with the formation of 
fog. The results for a +/-1 **C change in in the no 
fog cases, are a -/+4.3X variation in the moisture flux. 

Bea surface t*<np*rature was the next variable tested. 
The range in this case was from S to 14 °C. Sea surface 
temperature effects the heat flux directly and the 
moisture flux through q^. Figures 4.3.3 and 4.3.4 show 
the variations of heat and moisture fluxes, respectively. 

As shown in the figures, both the heat and moisture 
fluxes increase with increasing 0 q. For a +/-1 "C change 
of 0 q, the heat flux varies +/-3X, and the moisture flux 




Figure 4.3. 1 . 


Variation o-f Haat Flux with Cloud Top 
Tamparatura 
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v«rl«i +/-7X. 

Nsw t , tha intagratad liquid watar content was taatad 
ovar tha ranga of 300 to 480 m g/kg, for thia caaa. Tha 
intagratad liquid watar contant af facta tha haat flux 
through 0^ and tha moiatura flux through q^. Flgura 

4.3.5 ■howi tha variation of tha haat flux, and Flgura 

4.3.6 shows tha variation of tha moiatura flux. 

As can ba a man in Flgura 4.3.S, tha haat flux 
incraaaaa with incraaaing q^. Analysis shows that for a 
+/-I00 m g/kg changa in q 14 , tha haat flux changaa 
+/-0.25X. However, in Flgura 4.3.6, one may saa that tha 
moiatura flux dacraaaaa with incraaaing Q^. A +/-100m 
g/kg changa in q^ producaa a -/+3.3X changa in tha 
moiatura flux. 

Figuraa 4.3.7 and 4.3.8 ahow tha variation of tha 
haat and moiatura fluxaa with varying intagratad watar 
vapor contant. Tha intagratad watar vapor waa variad from 
720 to 900 m g/kg, for thia caaa. Tha haat and moiatura 
fluxaa ara affactad by q yi through 0j and q t , 
raapactlvaly. 

As shown in tha figuraa, tha haat flux dacraasas and 
tha moiatura flux incraasas with incraasing q yi . Analysis 
shows that for a ♦/-100m g/kg changa in Q yl , tha haat flux 
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Figure 4.3.5. Variation of Haat Flux with 
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Figure 4.3.6 


Variation of Moisture Flux with Integrated 
Liquid Water Content 
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Figure 4.3.7. Variation of Haat Flux with Intagratad Mata 
Vapor Contant 
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varin -/+0.6X and the moisture flux varies +/-SX. 

4.4 EXPECTED MODEL ACCURACIES WITH SATELLITE DATA 

The expected model accuracies with current satellite 
technology is now presented. Through infrared 

measurements f the cloud top temperature may be found 
within +/-2 "C. The accuracies of the sea surface 

temperature, the 10m neutral wind speed, and the friction 
velocity are +/-1.5 **C, +/-2m/n or 20X, and +/-20X, 
respectively (gfifcallitC Hiccowave Remote SffQsiQfl, edited 
by T.D. Allan, 19B3). These values were obtained from 
results of 8EA8AT and NIMBUS 7 data analysis. The 
expected retrieval errors of the integrated liquid and 
water vapor contents for the 8MMR package on NIMBUS S are 
+/-650m g/kg and +/-1500m g/kg, respectively (Chang and 
Wilheit, 1979). 

By applying these accuracies to the percent 

variations of the model discussed in the previous section', 

l 

the expected accuracy of the model may be found. Since 

the difference of the variations with respect to u . - and 

xun 

u # is IX, the larger of the two errors will be used. 

Qiven a 20X error in satellite derived u.^_ , the 

ion 

model will produce a 20X error in both the heat and 
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moiatura Huh—. A 2 °C trror in 0 ^ producaa a BX arror 
in tha haat flux and a 9X arror in tha moiatura flux. An 
arror of 1.3 °C in tha aaa aurfaca tamparatura givaa a 
4.5X arror in tha haat flux and an 6.75X arror in tha 
moiatura flux. A 630m g/kg arror in q^ producaa a 1.6X 
arror in tha haat flux and a 22.73X arror in tha moiatura 
flux. An arror of 1300m g/kg for q yi givaa a 9X arror in 
tha haat flux and a 73X arror in tha moiatura flux. 
Taking thaaa arrora to ba indapandant of aach othar, it ia 
poaaibla to computa tha axpactad r.m.a. arror of tha 
modal. For tha arrora liated abova for aach paramatar, 
tha axpactad ovarall r.m.a. arror of tha modal ia 24X for 
tha haat flux and 82X for tha moiatura flux (Saa Tabla 
4.4.1). Although thaaa arrora ara quita high, it muat ba 
undaratood that tha aatallita tachnology which allows tha 
maaauramanta of tha nacaaaary paramatara la atlll in ita 
infancy. Tha major arrora in tha modal may ba attributad 
to tha calculationa of q^ and q yi . Futura improvamant in 
tha algorlthma uaad to find thaaa paramatara would 
aignif leant ly improva thia modal *a aatlmataa of tha 


aurfaca fluxaa 



Table 

4.4. 

1. Individual Expected 
Satellite Parameter 

Model Error 

Errors due 
and Total 

to Each 
Expected 

PARAMETER 

8ATELLITE RMS ERRORS 

HEAT 

MOISTURE 

u 10n 


20X 

20X 

20X 

©b 


2.0 ~C 

8X 

9X 

o 

<D 


' 

U 

0 

in 

a 

H 

4.5X 

6.7SX 

*li 


650 m g/kg 

1.6X 

22.75X 

q vi 


1500 m g/kg 

9X 

75X 


Total 

Expected R.M.8. Errors* 

24X 

82X 




61 


4.3 STRENGTHS AND WEAKNESSES 

The modal '• strengths and weaknesses art now 
discussed. Tha modal handlaa wall in cloudy and cloudy 
with fog conditiona. Tha modal la extremely quick, and, 
given -reasonable input accuraciaa, will produca raaaonabla 
aatimataa of tha aurfaca fluxaa. Aa far aa tha author can 
tall, tha modal ia tha only ona currantly aval labia which 
will provida aatimataa of tha aurfaca fluxaa for a cloud 
toppad MABL from aatallita maaaurabla paramatara. 

On tha othar hand, tha modal ia not capabla of 
calculating tha aurfaca fluxaa whan a cloud ia not 
praaant. Tha calculation of tha fluxaa ia alao unclear 
whan gapa in tha clouda ara praaant. If tha gapa ara 
larga enough, tha calculation of tha cloud top temperature 
becomea difficult, and tha profilaa assumed in tha model 
break down. Tha modal 'a assumption that tha cloud top is 
also tha base of tha inversion ia not always correct. Tha 
cloud may extend into tha inversion layeri however, due to 
the uncertainty of how far into tha cloud tha infrared 
cloud top temperature ia actually being measured, 
this waaknasa may be corrected. Tha presence of fog below 
tha cloud layer also presents a problem due to tha 
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uncartainty of how tha modal 'a aquations raact to fog. 
Tha main waaknaaa of tha modal ia tha paramatarization of 
tha watar vapor abova tha boundary layar. Tha maan valuaa 
chosan ara not rapraaantati va of all poaaibla caaaa, and 
furthar work in this araa ia naadad. 



CHAPTER 3 i CONCLUSIONS 


In conclusion, this thssis has produced a viable 
model to estimate the surface fluxes from satellite 
measurable parameters. The model applies to cases of cold 
air outbreak, California coastal stratus, or to any cloud 
topped MABL which has a sufficiently strong heat flux or 
cloud top radiative cooling to cause the potential 
temperature and total water profiles of the layer to be 
well -mixed. With improvements of the current satellite 
technology, this model will provide accurate estimates of 
the surface fluxes. The model is relatively simple and 
extremely quick running. The model is also the only one 
of its kind currently available to estimate the surface 
fluxes of momemtum, heat, and moisture of the cloud topped 
marine atmospheric boundary layer from satellite 
measurable parameters. 

In light of this thesis, some areas of further 
research are suggested. More work is necessary to find a 
method of parameterizing the water vapor content above the 
MABL. More research is also necessary for testing the 
model for the case of California coastal stratus to see 
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how wel 1 it reacts under nearly neutral conditions. The 
evaluation of the surface fluxes under foggy conditions 
and how well the model handles them is another area for 
study. The problem of gaps in the cloud deck and a method 
to handle this is an area for further thought. Also, 
other possible sets of parameters should be examined to 
see if a better method than the one presented may be 


found 
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